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We use the adaptive time-dependent density matrix renormalization group method (t-DMRG) to 
study the nonequilibrium dynamics of a benchmark quantum impurity system which has a time- 
dependent Hamiltonian. This model is a resonant-level model, obtained by a mapping from a 
certain ohmic spin-boson model describing the dissipative Landau-Zener transition. We map the 
resonant-level model onto a Wilson chain, then calculate the time-dependent occupation nd(t) of 
the resonant level . We compare t-DMRG results with exact results at zero temperature and find 
very good agreement. We also give a physical interpretation of the numerical results. 

PACS numbers: 03.65.Yz, 74.50.+r, 33.80.Be, 73.21. La 
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I. INTRODUCTION 

Quantum impurity models, describing a discrete de- 
gree of freedom coupled to a continuous bath of exci- 
tations, arise in many different contexts in condensed 
matter physics. In particular, they are relevant for the 
description of transport through quantum dots and of 
qubits coupled to a dissipative environmemi 2 .. In recent 
years, there has been increasing interest in studying the 
real-time dynamics of such models for Hamiltonians H(t) 
that are explicitly time-dependent, as relevant, for exam- 
ple, to describe external manipulations being performed 
on a qubit. It is thus important to develop reliable nu- 
merical tools that are able to deal with such problems 
under very general conditions. 

The most widely used numerical method to study 
quantum impurity systems is Wilson's numerical renor- 
malization group (NRG)^. With the recently proposed 
time-dependent NRG (TD-NRG)^ one can now calculate 
certain class of time-dependent problems where a sudden 
perturbation is applied to the impurity at time t = 0. 
TD-NRG may very well be accurate for arbitrary long 
time. However, up to now, TD-NRG is not capable of 
dealing with a Hamiltonian H(t) with a time-dependence 
more general than a single abrupt change in model pa- 
rameters at t — 0. We will show in this paper that the 
adaptive time-dependent density matrix renormalization 
group method (t-DMRG) is a promising candidate for 
treating a general time-dependent Hamiltonian Hit). 

The density matrix renormalization group method 
(DMRG) is traditionally a numerical method to study the 
low lying states of one-dimensional quantum systems 5 . 
The recent extension of this method, the adaptive 
time-dependent DMRG (t-DMRG) 6 ^, can simulate real- 
time dynamics of one-dimensional models with time- 
dependent Hamiltonians as well. t-DMRG has already 
been used to study problems involving real-time dy- 
namics of one-dimensional quantum systems, for exam- 



ple the far-from-equilibrium states in spin- 1/2 chains 8 , 
dynamics of ultracold bosons in an optical lattice*^, 
transport through quantum dots^i, dynamics of quan- 
tum phase transition 12 , and demonstration of spin charge 
separation^ 3 .. These works showed that t-DMRG is a 
versatile and powerful method to study the real-time dy- 
namics of one-dimensional quantum systems. 

The underlying mathematical structures of DMRG and 
NRG are similar in the matrix product state represen- 
tation language-^. Indeed, once a quantum impurity 
model has been transformed into the form of a Wilson 
chain model, it can be treated by DMRG instead of 
NRG 1 - ' 15 ^ 6 ' 17 . This possibility opens the door toward 
studying time-dependent quantum impurity models us- 
ing t-DMRG. In this paper, we take a first step in this 
direction by using t-DMRG to study a simple, exactly 
solvable quantum impurity model whose Hamiltonian is 
a function of time. This model allows us to benchmark 
the performance of t-DMRG by comparing its results to 
those of the exact solution. 



II. THE MODEL AND DMRG METHOD 

We study the resonant-level model with a time- 
dependent potential applied to the level. The Hamil- 
tonian is 

H(t) = e d (t)dU + J2£kc{c k + Vj2(d f c k +c{d). (1) 



fff creates a spinless fermion on the level (impurity) and 
ct creates a spinless fermion with momentum k in a con- 
duction band whose density of states is constant between 
—D and D and zero otherwise, with Fermi energy set 
equal to 0. The energy of the local band is swept linearly 
with time, ed(t) = Dvt, where v is the sweeping rate in 
units of the half band width D. This model is equiva- 
lent to the dissipative Landau-Zener model with a Ohmic 
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boson bath whose spectral function is J(u>) = 27raw, for 
lu <C u) c , where lo c is the high energy cutoff^, and the di- 
mensionless strength of dissipation parameter a is hence- 
forth set equal to \. When a is close but not equal to 
i, the Hamiltonian (1) contains an additional interaction 

term proportional to U(d^d — ^)(J2k v c k c k' ~ \)r^i but 
this case will not be considered here. 

At time to — » — oo the local level contains a spin- 
less fermion and the band is half filled. Then, we lift 
the energy of the level linearly with time. As the level 
approaches the band, the probability that the fermion 
jumps to and from the band will increase, and decrease 
after the level has passed the band. In this paper we 
study this problem in detail. In particular, we are inter- 
ested in the expectation value of the occupation number 
on the level nd(t) at time t. 

Before using t-DMRG to solve this problem, we need 
to transform the Hamiltonian to a DMRG-friendly form. 
This can be realized by using a standard Wilson map- 
ping (originally invented in the context of NRG), which 
include two steps: logarithmic discretization of the band 
and converting the Hamiltonian to a hopping for m 20 ' 21 . 
Here, we just give the final result: The Hamiltonian ([I]) 
is mapped to a semi-infinite Wilson chain 

H(t) = e d (t)dU + ( — + dVe>) 

7T 

oo 

+ -(1 + A" 1 ) A-^Uftfn+i + fl+ifn), (2) 

wherein = (1 - A _n_1 )(l - A -2n-1 )~5 (l - A~ 2 "~ 3 )~5. 
r = npV 2 is the hybridization parameter, and p is the 
density of states at the Fermi level. A > 1 is a log- 
arithmic discretization parameter, which means we di- 
vide the band into discrete energy intervals determined 
by ±A _1 , ±A~ 2 , ±A~ 3 , • • ■ . In the limit A -> 1, the dis- 
cretized spectrum becomes dense throughout the band. 
The hopping factors in Hamiltonian @ decrease expo- 
nentially, so it is sufficient to keep the first L sites to 
achieve an energy resolution of A~ L / 2 . 

The dimensionless parameter r = 2T/v can be used to 
define three typical regimes of this problem. They are: 

• Fast sweep: r«l 

• Intermediate sweep: r ~ 1 

• Slow sweep: r» 1 

We will examine the performance of DMRG in all these 
regimes. 

The Wilson-chain form of Hamiltonian ([2]) can now 
be treated using DMRG. We first use infinite and fi- 
nite DMRG^ to calculate the ground state of the initial 
Hamiltonian H(to) at to. This ground state is a very good 
approximation to the true initial state in the ideal case 
in which the level would start from to — * — oo as long 
as e<i(to) <C — |r|. In the fast and intermediate sweep 
regimes, we can choose to so that the e<i{to) — Dvto is 



far below the Fermi surface to satisfy e d (to) *C — |T|. In 
slow sweep regime we can do the same if we use a very 
large |to|- However, a more efficient way we adopt is to 
use a moderate toj but set ed(to) as a very low value (e.g. 
— 10000-D). After we get the starting state we apply the 

evolution operator Te 1 H ^ ds on the starting state 
|*(f )) to get the state |*(t)) at time t using t-DMRG: 

=Te- l/ <o " {s)ds \y(t Q )). (3) 

Here T is the time-ordering operator, and we set h = 1 
in this paper. 

More specifically, we first divide the time interval t 
into a series of tiny time steps of the length r. The 
Hamiltonian is a function of time, but in each tiny time 
step it can be approximated by a constant, so we have 

7- e - ,; /4 "^ ds ~ e -*r£(t-5) . . . e -trJff(|T) e -trd(5) < ( 4 ) 

We chose the the value of Hamiltonian in the middle of 
each interval to represent the Hamiltonian of that inter- 
val. At every time step we decompose e - lH ( s ) T into local 
operators using second order Suzuki- Trotter decomposi- 
tion, and we get 

e -iH(s)r _ e -iT[H d . (s) + HoA+H la + --- + H L . liL ] 
_ g— if -H"<i,o(s)g— if Ho.ig— if -Hi, 2 
. g — Hl<— i,L g— if Hi* — 1,L 
e-if#i,2g-ifHo,i e -ifi? <i ,o + (t 3 ), (5) 

where 

H dt0 (s) =e d (s)dU+(—)i(f o i d + rff ), (6) 

and H n/n+ i is the hopping term involving site n and n+1. 
The only time-dependent part of the Hamiltonian is the 
impurity, so we only need to update the Suzuki- Trotter 
term of the impurity and the first site of the Wilson chain 
e -2ffld,oO) a t every time step. 

We can also easily extend this method to study fi- 
nite temperature dynamics. Instead of using infinite and 
finite DMRG to find the starting state, we use finite- 
temperature DMRG 2 ^ to get the starting state. Then, 
one can evolve this purified state using t-DMRG to sim- 
ulate the real-time dynamics at finite temperature^ 3 -. In 
this paper, however, we only focus on the zero tempera- 
ture and noninteracting case. 

III. EXACT METHOD 

The Hamiltonian @ is of quadratic form, so we can 
write it as 

H(t) = (at, a|, • • • , a^_ 1 ) J ff(t)(a , a u ■ ■ ■ , a L _i) T , (7) 

where a = d, a.j = H{t) is a L x L Hermitian 

matrix with L being the length of the Wilson chain. 
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By diagonalizing H(to) we get 

H(t ) = y^E k ala k . 
k 

The fc-th single particle state is 

ifc) = 4i°>=E^ a li°)' 



(8) 



(9) 



where Ui k are the eigenvectors of i/(to), in the sense that 
J2j H(to)ij u jk — E k u ik . 

At to the system is in its ground state, characterized 
by the single-particle distribution function 

«*>-{!: S<S ■ (10) 

The initial density matrix of the whole system is 

p(t )-E/( fc )l fc )( fc l- ( U ) 

k 

The density matrix evolves according to the von Neu- 
mann equation 



(12) 



This equation can easily be solved with an ordinary dif- 
ferential equation solver such as Matlab's ode45. Then 
we can calculate the expectation value of operators, like 
hd(t), as 

n d (t) = Tr [n d p{t)\ = Tr aja yo(i) . (13) 



IV. RESULTS AND PHYSICAL 
INTERPRETATION 

In Fig. [1] we plot both the exact and DMRG results in 
the three typical parameter regimes at zero temperature 
respectively. We use Wilson-chain length L = 160 and 
logarithmic discretization parameter A = 1.08 for all the 
three figures. We will discuss the discretization method 
in more detail in the next section. Note that we set D = 1 
in our calculation. 

For all three regimes, the DMRG error (shown in Fig. 
[2]for fast regime) is at worst of order 10~ 4 when keeping 
X = 100 states during DMRG calculation. This error can 
be further reduced by increasing \. 

Let us now try to understand the results physically. In 
the fast sweep regime the spinless fermion on the impu- 
rity does not have enough time to totally jump into the 
band, so the occupation on the impurity rid(t) converges 
to a finite value as the level is swept through and out 
of the band. In contrast, in the slow sweep regime the 
fermion ends up in the band with a very high probability. 
For comparison we also show the results of an adiabatic 
sweep in the slow sweep regime in Fig. [TJ The adiabatic 
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FIG. 1: (Color online) The local occupation number nd(t) as 
a function of time, calculated with both exact and DMRG 
method in the three parameter regimes. At the top, we give 
the choices made for the following parameters: sweeping speed 
v, Suzuki- Trotter step r, logarithmic discretization parameter 
A, Wilson-chain length L, and the number of states kept in 
DMRG calculation \- The value of hybridization parameter V 
and the corresponding dimensionless parameter r = 2T/v are 
given in each figure respectively. The insets zoom in on fine 
details of the curves, (a) The markers in the inset indicate 
the periods of the oscillations and beats obtained from the 
simple physical picture discussed in the text [see Eq. 1)190 ]. 
(c) The dashed lines are the reference results of the adiabatic 
sweep calculation, and the dash-dotted line is the adiabatic 
sweep result in infinite band limit, which is simply rid(^) = 
i — — arctan ^ . 
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r=0.05D, v=lD, T=0.05/D, A=1.08, L=160 
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FIG. 2: (Color online) Error of the DMRG results for Sn d (t) = 
n™ RG (i) - n^ xact (£) in fast sweep regime when keeping 100, 
150 and 200 states. 



results are obtained from the thermodynamic average 
Tr[p ed / t ^nd], where p €d (t) is calculated using Eq. (fTTj) with 
single particle states \k ed ( t )) of the Hamiltonian H ed ( t y 
Evidently, the DMRG and exact results agree very well 
with the adiabatic results. 

Another important feature of the results is the oscilla- 
tion of nd{t). To understand it, we first study a simplified 
model, in which we only consider one level in the band 
and disregard the rest levels for the moment. When there 
is one spinless fermion in this system the Hamiltonian is 



Hit) 



E (t) 7 
7 Ex 



(14) 



This is just the Hamiltonian of the original Landau-Zencr 
problem. We denote the instantaneous two eigenstates as 
\+)t, \—)t with the corresponding eigenenergies E±(t) = 
\{E {t) +E X ±w(t)], where 



w(t) = V47 2 + (Ei -Eo(*)) 2 



(15) 



The probability that a state of the form \4>(t)) = 
a\—)t + b\+) t at time t will still be found in the same 
state at time t + St, is given by 

P{t) = Mt)\<f>(t + St))\ 2 (16a) 
= |a| 4 + \b\ 4 + 2\ab\ 2 cos [uj(t)St]. (16b) 

In each time interval, the instantaneous oscillation fre- 
quency u;(t) of P(t) is equal to the instantaneous oscil- 
lation frequency of |(</>(io)|</>(i))| 2 to the zeroth order in 
St. Therefore, the probability for the system initially in 
a state \4>(to)) to still be found in this state at a later 
time t, 



p(i) = |<<^oMt)}| 2 



(17) 



will have an oscillating component proportional to 

J t 



cos 



Stn L °^ dS 



We now return to the original problem and use the 
picture described above to roughly estimate the period 
of the oscillations in the fast sweep regime. In the fast 
sweep regime according to Pauli exclusion principle the 
influence of the unoccupied levels of the upper half of 
the band is dominant. We can neglect the lower half of 
the band, and add up the contributions of all levels E\ 
in the upper half band to the oscillations by integrating 
the above mentioned cosine term over the energies E\. 
Therefore the occupation on the resonant-level 



n d {t) 



P{t)dE-i 



(18) 



will contain an oscillating contribution proportional to 



~ D , 




jit -to) 


cos 



-(vt*-Dt + Dto - 



vti 



(19) 



To get the above result, we approximated u>(t) by E\ —vt, 
neglecting the term 4j 2 = 4TD/ir (7 = y/DT/n is the 
prefactor of the hybridization term in the energy repre- 
sentation of Hamiltonian ([I} 2 — •)■ This is a good approxi- 
mation, except around t — 0, when the local level is near 
the middle of the band, and \E — Ei \ is not significantly 
larger than T. 

The resulting Eq. (|19p can be used to understand the 
nature of the oscillations and beats observed in the fast 
sweep regime in Fig. [TJ The factor sin [D(t — to)/ 2] is 
the beat, and the period of the beats is Treats = 47r/D. 
We plot two markers with a separation of 4tt/D under 
the curve in the inset of Fig. [TJa); they fit the period of 
the beats very well. The markers above the curve in the 
insets of Fig. QTa) are obtained by solving 

i(ui 2 - Dt + Dt - vtl) = 2tott + const, (20) 

where m is an integer such that the markers are best 
aligned with the maxima of the oscillations shown. We 
can see that the final agreement in position is excellent. 

Last but not the least, we examined the dependence 
of the final local level occupation number n,j(+oo) on r 
(shown in Fig. [3]), and find it has the typical Landau- 
Zener exponential relation: 



n-di+oo) = e 



(21) 



This agrees with previous analytical result a 24 ' 25 . Note 
that though n^+oo) only depends on r, the detailed 
structure of the nd{t) curve is determined by v and T 
respectively. (See Eq. (|2"D|) for example.) 



V. ROLE OF DISCRETIZATION PARAMETER 

As in NRG, the value chosen for the discretization pa- 
rameter can affect the real-time dynamics, if it does not 
lie sufficiently close to 1. Fig. 0|a) compares the exact 
results of A = 1.08 and A = 2 in fast sweep regime. Note 
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TABLE I: The n d (+oo) data used in Fig. |3] 



r \ 


0.1D 


0.3D 


0.9D 


2.7D 


8. ID 


24.3D 


0.05D 


0.36516 


0.71499 


0.98419 


0.96345 


0.98766 


0.99587 


0.2D 


0.01831 


0.26199 


0.63925 


0.86164 


0.95155 


0.98358 


0.8D 




0.00480 


0.16758 


0.55114 


0.81984 


0.93593 


3. 2D 






0.00081 


0.09221 


0.45177 


0.76732 


6.4D 








0.00850 


0.20404 


0.58877 


12.8D 








0.00009 


0.04162 


0.34660 
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FIG. 3: (Color online) Exact results checking the rela- 
tion between the final local occupation number rid(+oc) 
and r. Numerically, we approximate nd(+oo) by averag- 
ing Ud(t) of the last 4 time steps. The time span we 
use here is t £ [-200/D, 200/D]. To get n d (+oo) at dif- 
ferent r, we choose 6 different F from a wide parame- 
ter regime, and with each F 6 different sweeping speed: 
v = 0.LD, 0.3D, 0.9-D, 2.7D, S.1D, 24.3D are used to calculate 
nd(+oo). We only plot the data for r < 10 because the ac- 
cumulated numerical error becomes significant compared to 
nd(+oo) for r > 10. The dashed line is a reference line of 
ln[nd(+oo)] = — r. The inset zooms in on small r. 
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that for A = 2, big oscillations in n<j(i) remain long af- 
ter the transition. These are artificial consequences of 
the rather coarse discretization scheme, which diminish 
strongly as A is reduced towards 1. Indeed, for A = 1.08, 
most of these oscillations have disappeared. Further re- 
duction of A does not change the results significantly any- 
more. Note that, incidentally, the ability of allowing a 
logarithmic discretization parameter very close to 1 is a 
big advantage of DMRG over NRG. 

With the physical picture described in the last section, 
we can also understand why there are artificial oscilla- 
tions if A is big. If we use a big logarithmic discretization 
parameter, the part of the band far away from the Fermi 
level is poorly represented by only a few levels, which 
means that the oscillations from different levels do not 
average out as well as would have been the case for a 



FIG. 4: (Color online) Comparison of the exact results of 
different logarithmic discretization parameters. Both figures 
zoom in on fine details. Here we study fast sweep regime 
as an example, (a) Comparison of the results with differ- 
ent A. The large oscillations for A = 2 for later times are 
due to the relatively coarse discretization and hence artifi- 
cial. (b),(c) Comparison of the converging speed with respect 
to the Wilson-chain length L of linear and logarithmic dis- 
cretization method. 

true continuum of levels. 

We use logarithmic discretization instead of linear dis- 
cretization because in the problem we studied, the levels 
near Fermi surface contribute more than levels far away 
from it, and logarithmic discretization represents the part 
of band around Fermi surface more efficiently,? 7 '. This is 
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reflected in the convergence of the results with respect 
to the Wilson-chain length L shown in Fig. 0J As other 
parameters are the same, the two discretization meth- 
ods will both converge to the same result when L — > oo. 
Therefore the faster the result converges the better the 
method is. We can see from Fig.|4](b) that the difference 
of Udif) between L = 120 and L = 160 chains is already 
negligible for the case of logarithmic discretization while 
still significant if using linear discretization, which means 
the results converge more quickly if we use logarithmic 
discretization. This is even more obvious by comparing 
the convergence speed of n^+oo) shown in Fig. SJc) . 

VI. CONCLUSIONS AND OUTLOOK 

By studying a benchmark model we demonstrated that 
the t-DMRG is a very accurate method to calculate real- 
time dynamics of quantum impurity system with a time- 
dependent Hamiltonian. To compare with the exact re- 
sults, the model we studied here is a non-interacting 
model, but DMRG can also treat interacting problems 
similarly. 

Though t-DMRG cannot calculate arbitrary long times 



(in contrast to TD-NRG) it can give reliable results in a 
relatively long time which we expect to be long enough 
for numerous practical purposes. For example, in quan- 
tum information, where fast quantum processes are more 
useful, the relevant physics happens in a relatively short 
time scale, which can be simulated by t-DMRG with a 
high precision. We thus expect t-DMRG to be a powerful 
tool to study the real-time dynamics of quantum impu- 
rity systems, in particular in the context of modeling the 
dynamics of damped, driven qubits. 
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